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We have analyzed the effects of the anisotropic energy bands of phosphorene on magnetoroton 
branches for electrons and holes in the two Landau levels close to the band edges. We have found 
that the fractional quantum Hall effect gap in the lowest (highest) Landau level in conduction 
(valence) band is slightly larger than that for conventional semiconductor systems and therefore 
experimentally observable. We also found that the magnetoroton mode for both electrons and 
holes consists of two branches with two minima due to the anisotropy. Additionally, we show that 
due to the anisotropy, there is a second mode with positive dispersion, well separated from the 
magnetoroton mode for small wave vectors. These novel features of the collective mode can be 
observed in resonant inelastic light scattering experiments. 


With the advent of graphene Hi the two dimensional 
(2D) materials have gained renewed attention during the 
last decade because of their remarkable electronic, optical 
and mechanical properties and for potential device appli¬ 
cations. Practical application of graphene in electronic 
devices is however limited due to its zero band gap. Al¬ 
though there are external means to create band gap in 
graphene, such as application of a bias voltage to bilayer 
graphene, in the past few years intense search has been 
undertaken to find other 2D materials, which in addition 
to exhibiting mobilities close to graphene, also posses a 
considerable band gap in their normal state. There were 
several Dirac materials, such as silicene Q, germanene 
[3|, M0S2 and other group-VI transition-metal dichalco- 
genides @ that have been lately explored and promising 
results were observed. 

In the past year another material which has gained 
popularity in this context, is the 2D version of black 
phosphorus (BP) the single layer of BP known 

as phosphorene. BP is the most stable allotrope of Phos¬ 
phorus at room temperature and pressure. Few-layer BP 
can be obtained through mechanical exfoliation method 
akin to graphene. However, unlike graphene the BP lay¬ 
ers are not perfectly flat. Due to the sp^ hybridization of 
3s and 3p atomic orbitals the layers of BP form a puck¬ 
ered surface. The band gap of BP depends on the number 
of layers of the sample, and it ranges from 0.3 eV in the 
bulk case at the Z point in the Brillouin zone to 1.5 - 2 
eV for few layer and monolayer case at the P point. The 
mobilities of the few-layer BP obtained so far lie in the 
range of 1000 - 5000 cm^V“^s“^ d,[HI, Ell- Considerable 
progress has been made in the study of the electronic and 
optic al prop erties of this material both experimentally 
[Q-S, lTll4T^ and theoretically [T3 - l^ . Higher mobilities 
were achieved by sandwiching the BP in hexagonal boron 
nitride flakes and placing on top of the graphite back 
gate that has enabled observation of the integer quantum 
Hall effect [l^. Clearly the next step in the experimen¬ 
tal progress is to reach the regime where the fractional 
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quantum Hall effect (FQHE) [2^ can be observed. 

The FQHE in Dirac materials shows many interesting 
features not found in conventional semiconductor quan¬ 
tum heterostructures. The FQHE states in monolayer 
and bilayer graphene were investigated theoretically [26l - 
[Hj | and experimentally Hili. It was shown that while 
the FQHE states in graphene in the n = 0 Landau level 
(LL) are completely identical to the case of conventional 
semiconductor systems, graphene exhibits more robust 
FQHE states in the n = 1 LL with a bigger gap than in 
the n = 0 LL (2^. That behavior is not present in con¬ 
ventional semiconductor systems. In bilayer graphene 
the application of a bias voltage results in some LLs 
the observation of the phase transition between several 
incompressible FQHE and compressible phases [2^, [Sll ■ 
The FQHE in silicene and germanene indicated that be¬ 
cause of the strong spin-orbit interaction present in these 
materials as compared to graphene, the electron-electron 
interaction and the FQHE gap are significantly modi¬ 
fied [13. The puckered structure of phosphorene ex¬ 
hibits a lower symmetry than graphene. This results in 
anisotropic energy spectra and other physical character¬ 
istics of phosphorene, both in momentum and real space 
in the 2D plane [^, [l^. Therefore it is desirable to ex¬ 
plore how this anisotropy of phosphorene manifests itself 
in the FQHE states. 

In this work we consider the FQHE states in phospho¬ 
rene for filling factor 1 / = 1/3 in the conduction band 
and ly = 2/3 in the valence band in two LLs closest to 
the band edges. We analyze the energy spectra and the 
magnetoroton modes of the system with finite number 
of electrons using the exact diagonalization scheme. We 
observe the FQHE gap in the lowest (highest) LL in the 
conduction (valence) band, with the gap being slightly 
higher than in conventional semiconductor systems. This 
indicates that the FQHE should be observable in phos¬ 
phorene. We also show that both the gap of the FQHE 
and the magnetoroton mode are similar for the lowest 
LL in the conduction band and highest LL in the va¬ 
lence band. The anisotropy of phosphorene manifests 
itself in two different ways. We find that the magne¬ 
toroton mode in phosphorene (for electrons and holes) 
is split into two branches with two minima due to the 
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anisotropy. We also find that the anisotropy causes a sec¬ 
ond mode with positive dispersion to appear that is well 
separated from the magnetoroton mode for small wave 
vectors and is therefore experimentally more accessible 
than for isotropic systems. For the FQHE states in the 
second LLs of the conduction and valence band of phos- 
phorene, the FQHE gap is an order of magnitude smaller 
than for the lowest LLs and approaches a soft mode, with 
possible incompressible to compressible phase transition. 

In what follows we consider a finite-size system of elec¬ 
trons in 2D phosphorene in the toroid al g eometry us¬ 
ing the exact diagonalization procedure [25| . The many- 
body Hamiltonian of the system is 


n = + (1) 

i i^j 


where "Hp is the single-particle Hamiltonian in phospho¬ 
rene and the second term is the Coulomb interaction. 
Since we are interested in the FQHE in the low lying 
LL, which are situated close to the conduction and the 
valence band edges, we employ the continuum approx¬ 
imation of the single-particle Hamiltonian derived from 
the microscopic tight-binding approach M- In all our 
calculations we take the x axis along the armchair direc¬ 
tion and the y axis along the zigzag direction. Then the 
single-particle Hamiltonian in single-layer phosphorene 
in the continuum approximation has the form [l9| 


= ( Hf(k) ) • < 2 ) 

Hoik) = Uq + + Vyk^, 

iLi(k) = (5 -f + -fykl + ixk^, 


where k is the wave vector and all parameters in Llo(k) 
and (k) are related to five hopping amplitudes in 
the tight-binding model [l^: Uq = —420meV, = 

10.1 meV • nm^, r]y = 5.8 meV • nm^, 6 = 760 meV, Xx = 


38.3 meV • nm^, jy = 39.3 meV • nm^, x = 525 meV • nm. 
For these parameters the band edges of the conduc¬ 
tion and the valence band are located at the energies 
= 340 meV and = — 1180meV, with energy gap 
Eg = 1520 meV. In a magnetic field, we first write 
k = —iV + eAjhc, where the vector potential is in the 
Landau gauge A = Bxy. For all the FQHE states con¬ 
sidered here the many-body ground state is fully spin po¬ 
larized. We henceforth disregard both the electron spin 
and the Zeeman energy. For the next step we cast the 
Hamiltonian @ into the second quantized form using 

y/^k^ - i^ky , where 


the ladder operator a = 


= “^y/VxVyl^o ^0 = y/kcJeB is the magnetic 
length. It is easy to verify that [a, a^] = 1. Then the 


components of the Hamiltonian ([2]) are 
Hq = uq + -I- , 

E[^ = 5 huj^ + 2 ^ ~ ~ 

-I- ikuJcX y 

where s = Xxhxy * = ly/Vy, X = xl‘^\/vj^c- We search 
for the eigenvalues of the Hamiltonian (I2|) in the form 




En=“ (x, y) \ 

YZ=t Bn(l)n,k^ix,y) j ’ 


( 3 ) 


where and Bn are constants and 4>n,k^ (x, y) are the 
usual harmonic oscillator wave functions 

oh 

<Pn,k^{x,y)= =e‘^yye 2 

^2^n\Lyy/l^ 

X Hn {y/a {x + ky£l)) , ( 4 ) 

where a = hu!^/2r]^, Ly is the sample size in the y direc¬ 
tion and ky is the wave vector in the y direction, which is 
conserved in the chosen gauge. Diagonalizing the Hamil¬ 
tonian matrix obtained with the wave function m we 
determine the eigenvalues of the Hamiltonian ([5]), which 
are the LLs of phosphorene and the appropriate wave 
functions by determining the constants A„ and Bn- 

In order to consider the many-body states we proceed 
in the same way as for conventional semiconductor sys¬ 
tems m, . We consider the many-body states in the 
system in a toroidal geometry, i.e., we apply the periodic 
boundary conditions (PBC) on both directions of the sys¬ 
tem m. In, m and cast the wave functions @ and (HD 
into the PBC preserving form. This process naturally 
defines the magnetic translation unit cell, which we take 
rectangular in our calculations with the sides and Ly, 
L^Ly = 27r£gAj,, where is the number of magnetic 
flux quanta passing through the unit cell and is an inte¬ 
ger number. also characterizes the degeneracy of each 
LL. We take = pN and = qN, where is the 
number of electrons, N is the greatest common divisor of 
Ag and N^, p and q are integers. The filling factor in a LL 
is defined as iz = p/q. We then construct the many-body 
states from the appropriate single-particle states and cal¬ 
culate the Hamiltonian matrix of the form © with this 
many-body basis states. We make use of the many-body 
translational symmetry to reduce the size of the Hilbert 
space We define the relative translation operator 

Tb(a) and note that the translations T^{pLjnn) preserve 
the Hilbert space and commute with the Hamiltonian ([!]) 
and with each other. Here L^n = "rnEx + nLy is the 
magnetic translations lattice vector, L^, and Ly are the 
unit vectors of magnetic translations unit cell, m and 
n are integers. Therefore we use these relative transla¬ 
tion operator (Tb(pLm„) eigenstates to bring the com¬ 
plete Hamiltonian matrix into the block-diagonal form 
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FIG. 1: (a) Magnetic field dependence of few low lying LLs in 
phosphorene, both in the conduction (blue) and the valence 
(red) band, (b) Magnetic field dependence of the v = XjZ 
gap in the conduction band and the ;/ = 2/3 gap in the va¬ 
lence band lowest LLs of phosphorene. The finite-size system 
considered for FQHE gap has six electrons. 


and also to characterize each of the many-body states 
with its appropriate relative momentum Icr. 

In Fig. [T] (a) the magnetic field dependence of the few 
low lying LLs both in the conduction band and in the 
valence band of phosphorene is shown. The dependence 
is almost linear and is similar to the case of conventional 
semiconductors. In fact as was shown previously by sev¬ 
eral authors [l^ [Tsl - f^ in the low energy regime the en¬ 
ergy bands of the conduction and the valence band can 
be approximated by parabolic dispersion with anisotropic 
masses in the x and y directions. Therefore for low-lying 
LLs and for magnetic fields up to ^25 T the dependence 
of the LLs on the magnetic field is linear. For this reason, 
in considering the FQHE states in the low-lying LLs we 
can disregard the LL coupling (the possibility for elec¬ 
trons to occupy states in higher LL). In case of systems 
with parabolic dispersion the LL mixing parameter At, 
the ratio between the Coulomb interaction energy and 
the cyclotron energy has the oc X/Vb dependence [s^. 
Therefore for high magnetic fields k < 1 and LL mixing 
can be safely ignored in our work. 

In Fig. [T] (b) the magnetic field dependence of the 
FQHE magnetoroton gap is presented for the lowest 
LL in the conduction and in highest LL in the valence 
band. Due to the particle-hole symmetry the filing fac¬ 
tor Az = 1/3 of electrons in the valence band corresponds 
to jz = 2/3 filling for holes. The results are presented 
for the finite-size system with = 6. The aspect ratio 
of the magnetic unit cell A = L^jLy is A = iVg/4 in all 
calculations of the FQHE states. As the gap is presented 
in units of = e^/e£g, where e = 10.2 is the dielectric 
constant, this dependence is due to the mixing of each LL 
in phosphorene in terms of the harmonic oscillator wave 
functions (|4|) . In conventional semiconductors the FQHE 
gap in units of does not depend on the magnetic field 
and for a similar system size the gap is 0.041. Therefore 
as seen in Fig. [T] (b) the FQHE gaps in phosphorene are 



FIG. 2: The energy spectrum of the six-electron FQHE sys¬ 
tem in phosphorene with = A in the lowest LL of the con¬ 
duction band. The red and blue dash lines show schematically 
the two branches of the magnetoroton mode. Green dash line 
shows the second mode described in the text. Inset: the mag¬ 
netoroton mode and the second mode for = 7. 

slightly larger than in conventional systems. The FQHE 
gap for the LL in the valence band is in fact larger than 
the gap in the conduction band, but overall all these gaps 
fall within the same energy range. The dependence of the 
FQHE gaps on the magnetic field is linear, albeit barely 
discernable, which again indicates that up to 25 T the 
non-trivial character of the LLs in phosphorene compared 
to conventional systems with parabolic dispersion is less 
important. The main features of phosphorene that can 
distinguish it from conventional semiconductor systems 
in the context of FQHE are the anisotropic nature of its 
bands, which we address next. 

Figure [2] shows the energy spectrum of the six-electron 
FQHE system in the lowest LL of the conduction band 
for phosphorene with filling factor u = 1/3. The inset 
shows the low energy collective modes for a seven elec¬ 
tron system. As in conventional semiconductor systems, 
the magnetoroton mode can be clearly identified. The 
effect of anisotropy of the energy bands of phosphorene 
is that the magnetoroton mode is split into two branches 
with two minima. These are schematically shown in the 
figure with red and blue dashed lines. These kind of 
branches were actually noted previously in the context 
of adding a band mass tensor (Galilean metric) to the 
Hamiltonian in order to investigate the correlation be¬ 
tween external and intrinsic metrics of the system 
[ 3 ^. The intrinsic metric was treated as a variational 
parameter which characterizes the shape of the correla¬ 
tion function of the Laughlin state and is obtained by 
finding the maximum overlap of the exact wave function 
and the family of Laughlin wave functions parametrized 
by the intrinsic metric. The splitting of the magnetoro¬ 
ton mode into two branches is possible in graphene in 
the case of anisotropic electron-electron interaction [s^. 
The appearance of two branches can be attributed to the 
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FIG. 3: The energy spectrum of the six-electron FQHE sys¬ 
tem with = I in the second LL of conduction band for 
phosphorene. Inset: the same spectrum but for a conven¬ 
tional semiconductor system (such as GaAs). 

fact that in the anisotropic system the minimum of the 
magnetoroton mode appears in different \kp\ for differ¬ 
ent direction of the relative momentum From this 
consideration we can conclude that phosphorene presents 
a real anisotropic system where the ideas of correlation 
of the external metric (mass tensor in this case) and the 
intrinsic metric can be explored. 

Interestingly, in Fig. we also observe a second mode 
(schematically shown by a green dashed line) for small 
wave vectors. It has positive dispersion and merges with 
the continuum for wave vectors > 1. This mode was 
observed previously in an experiment on the GaAs sam¬ 
ple and was attributed to a two-roton state [d^. One 
possible theoretical explanation for the appearance of 
the second mode was that two types of shear forces ex¬ 
erted on the Imuid state have different direction for dif¬ 
ferent modes [dy . While for isotropic systems these two 
modes should converge to the same point at —>■ 0 our 
studies indicate that for an anisotropic system these two 
modes are clearly separated for small wave vectors and 
do not converge to the same point as k^ —?> 0. There¬ 
fore phosphorene can greatly simplify the identification 
of the second mode from the magnetoroton mode in res¬ 
onant inelastic light scattering experiment. It should be 
noted that possibly anisotropy has played some role in 


the identification of the second mode even in the original 
experiment [dO ). In that experiment the applied mag¬ 
netic field is not perpendicular to the sample. The effect 
of the tilted magnetic field has been exhaustively stud¬ 
ied in the experiments on FQHE and it was shown that 
tilting results in anisotropic transport properties of the 
sample [d^. The effect of the tilted field on the FQHE 
states was also addressed theoretically by several authors 
and it was found that inclusion of the tilted field is 
similar to adding an anisotropic metric in the system, al¬ 
though for the tilted magnetic field the dynamics is more 
complicated than simply the inclusion of the anisotropic 
mass. Therefore, due to the tilted field the system per- 
haps acquire some sort of anisotropy in the experiment 
[4(]| . which reveals the second mode. Einally, we have 
performed similar studies for the highest LL in the va¬ 
lence band for = 2/3 and obtained very similar results. 

In Fig. [3l we present the energy spectrum of the six- 
electron FQHE system in the second LL of the conduc¬ 
tion band for phosphorene with filling factor u = 7/3. 
The same result for the conventional semiconductor sys¬ 
tem (such as GaAs) is shown as inset. It seems that the 
two branches of the magnetoroton mode and the second 
mode are still present. The gaps are however, very small 
and the system seems to be approaching a soft mode, re¬ 
sulting in a phase transition from the incompressible to 
a compressible phase. Similar results are also found for 
the second LL in the valence band. 

To summarize, we have studied the FQHE states in 
phosphorene in two LLs close to the band edges in the 
conduction and valence band. We have shown that the 
gaps for V = 1/3 of electrons in the lowest LL of the 
conduction band and = 2/3 for holes in the highest 
LL of the valence band are slightly larger than for con¬ 
ventional semiconductor systems. Therefore the FQHE 
states should certainly be observable in phosphorene. 
The anisotropic band structure of phosphorene causes 
splitting of the magnetoroton mode into two branches 
with two minima. Eor small wave vectors, there is also 
a second mode with positive dispersion that is clearly 
separated from the magnetoroton mode. These features 
should be observable in inelastic light scattering experi¬ 
ments. 
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